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Abstract 

A straightforward procedure has been developed to 
quickly determine an in viscid design of a hypersonic 
wind tunnel nozzle when the test gas is both calorically 
and thermally imperfect. This real gas procedure divides 
the nozzle into four distinct parts: subsonic, throat to 
conical, conical, and turning flow regions. The design 
process is greatly simplified by treating the imperfect gas 
effects only in the source flow region. This 
simplification can be justified for a large class of 
hypersonic wind tunnel nozzle design problems. The 
final nozzle design is obtained either by doing a classical 
boundary layer correction or by using this inviscid design 
as the starling point for a viscous design optimization 
based on computational fluid dynamics. An example of a 
real gas nozzle design is used to illustrate the method. 
The accuracy of the real gas design procedure is shown to 
compare favorably with an ideal gas design based on 
computed flow field solutions. 

Nomenclature 
Ci = constant in Eq. (I) 

C (i = contraction coefficient 

M = Mach number 

r* = nozzle throat radius or height 

R c = throat radius of curvature/r* 

5 = R ( + 1 

x = axial coordinate 

X a = start of subsonic contour with Eq.( 1 ) 

y = radial coordinate 

)\ = nozzle wall radial coordinate 

y M , = nozzle wall slope 

y = gamma, ratio of specific heats 

9 a = subsonic approach angle 

0 SU per = nozzle inflection angle 

Introduction 

A number of hypersonic wind tunnel facilities have 
required new wind tunnel nozzles either to replace 
existing nozzles or to extend the facility capabilities into 
new flight regimes. These changes have been driven in 
part by the effort to develop a viable scramjet engine for 
hypersonic vehicles and for access-to-space initiatives. 
To simulate this type of vehicle and engine environment 
at hypersonic velocities requires stagnation conditions 
where the test gas typically behaves both as calorically 
and thermally imperfect through part of or the entire 


nozzle. The inviscid design of a supersonic (and 
hypersonic) wind tunnel nozzle for a perfect gas was 
originated by Busseman in the 1920s (as noted in Ref. 1) 
and, in the opinion of the author, was finally perfected by 
Sivells 2 in the 1960s. No similar progress has been 
made for real gas wind tunnel nozzle design. This paper 
concentrates on an approximate inviscid nozzle design 
process for a calorically and thermally imperfect test gas. 

Supersonic Wind Tunnel Nozzle Desig n 
Ideal Gas Design 

The theory used in the design of supersonic wind 
tunnel nozzles for an ideal gas can be found in a number 
of reports 1 6 and books. 7 ' 9 The basic assumption made is 
that the boundary layer thickness is small compared to the 
characteristic length (nozzle radius), so the nozzle flow 
field can be treated as inviscid for designing the 
aerodynamic lines. Once the aerodynamic lines are 
determined, a correction is made to account for the 
displacement thickness of the boundary layer. This basic 
procedure has been applied successfully to many 
supersonic and hypersonic nozzles. A computational 
fluid dynamics (CFD) analysis should then be made to 
validate the flow quality and the design Mach number for 
the nozzle design. If the desired flow quality or Mach 
number is not obtained, this failure usually means that the 
boundary layer assumption was invalid and a different 
design method, which directly includes the viscous 
effects should be used. 10 



Fig. 1 Wind tunnel nozzle flow regions. 

The classical inviscid design procedure requires the 
position of the sonic line at the throat, a prescribed 
centerline Mach number, and the condition of uniform 
flow at the nozzle exit (Fig. 1). Sivells 2 determined the 
appropriate centerline Mach number distribution that 
would result in a nozzle wall with continuous wall slope 
and curvature from the throat to the exit. The centerline 
Mach number distribution for the upstream and turning 
sections is based on polynomial functions that match the 
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flow field solutions at each flow region boundary. This 
distribution results in a wall contour shape that is easier to 
machine. 

Once a centerline Mach number distribution has been 
specified, the classical design procedures use an inverse 
method-of-characteristics (MOC) design procedure to 
determine an inviscid contour. This inviscid contour can 
then be scaled to any particular size. The design contour 
for one operating condition is determined by computing 
the displacement thickness by using a boundary layer 
(BL) code and adding it to the inviscid contour. Modem 
implementations of MOC-BL design procedures typically 
differ in the way the sonic line and centerline Mach 
number distributions are specified and in the numerical 
techniques used to compute the inviscid flow field and the 
boundary layer solution. 

Real Gas Design 

Real gas effects are important in the design of many 
hypersonic blow down and reflected shock facilities. 
Recent advances in wind tunnel nozzle design with CFD 
codes coupled to optimization techniques" or with 
traditional MOC-BL design procedures are two options 
for designing a nozzle. A CFD-based method is the more 
accurate design procedure when the boundary layer is 
extremely thick. 10 However, for reflected shock tunnels, 
the run time is short enough that the nozzle wall stays 
cold, and the boundary layer is thin enough that accurate 
designs can be obtained with the classic MOC-BL 
approach. The CFD solutions for the reflected shock 
nozzles are extremely time-consuming due to the fine grid 
requirements for accurately simulating the thin nozzle 
boundary layers. 

At least two approaches have been used to include real 
gas effects in an MOC design procedure. The first 
approach is to include real gas behavior by using 
tabulated thermodynamic properties, assuming a 
thermally perfect gas. 12 The second approach is to divide 
the nozzle into flow- regions, which include a source flow 
region where the source flow is computed with a one- 
dimensional real gas code and is connected to a turning 
contour computed with an MOC design procedure for a 
calorically perfect or imperfect gas. 1314 Neither method 
addresses high pressure effects nor exploits a centerline 
Mach number distribution which guarantees a nozzle wall 
with continuous wall curvature. Both designs methods 
start with a source flow region at the nozzle throat. 

A variation of the second approach is used to develop a 
design procedure applicable to a real test gas with both 
caloric and thermal imperfections that would produce a 
wall design with continuous wall curvature. The method 
described in this paper differs from previous methods in 
that the approach follows Si veils' design philosophy of 
prescribing a centerline Mach number distribution and 
includes high pressure effects in the throat region. The 
author has successfully used this method to design a 
number of hypersonic wind tunnel nozzles at NASA 
Langley Research Center and elsewhere. 15 


Procedure for Real Gas Inviscid Hypersonic Nozzle 

Desig n 

The design procedure will be described assuming that 
the nozzle will have an axisymmetric geometry. The 
process and the computer codes developed can also 
handle two-dimensional nozzle design. The inviscid 
contour design is obtained by dividing the nozzle into 
four sections: (1) subsonic and throat section, (2) throat 
to conical section, (3) conical or source flow section, and 
(4) turning contour downstream of the conical section. 
Five specific steps are used to design the complete 
inviscid contour. For describing the design method, we 
will assume that the throat radius, throat radius of 
curvature ratio, stagnation conditions, compressibility 
factor, specific heats ratio (gamma), and gas composition 
are known. (Sometimes the exit diameter is given and the 
throat radius needs to be determined.) 

The first four steps will be completed with 
coordinates made nondimensional by the throat radius. 
The final design is obtained by scaling the nozzle design 
to obtain the desired throat radius or exit radius. 

The procedure will be explained and demonstrated 
with various options in Si veils’ design code 1 * for 
designing the inviscid contour of nozzle sections. In 
addition, an approach for defining a centerline Mach 
number distribution for use in a CFD-based design 
process will also be described. 

Once the inviscid design is completed, a BL code 
(like Anderson’s) 17 can be used to correct the nozzle for 
the displacement thickness, assuming turbulent flow and 
chemical and vibrational equilibrium. The nozzle design 
can then be validated with computed CFD solutions by 
using the procedure and the codes described by Korte. 18 
Typically the final nozzle contour length will be 
determined by truncating the nozzle design based on the 
boundary layer growth observed in the CFD flow field 
solutions or from facility length restrictions. 



Fig. 2 Subsonic and transonic nozzle geometry. 


Step 1 : Subsonic-Throat Region 

The geometry of a continuous curvature throat section 
used in the design method is given in Fig. 2. To connect 
the subsonic side of the throat section to an entrance cone 
or wedge, the subsonic contour was specified by the 
equation 1 9 
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>’„■ / r* = C, [l - exp(- (jr / r *) 2 /( 2C , R ))] + 1 ( I ) 

where 

C, = R { tan 2 (0Jexp(l) (2) 

and the point of intersection with the entrance cone or 
wedge is 

x„ / r* = - R c tan ( 6 a )-y/exp(l) . (3) 

This subsonic contour has the desirable qualities of 
providing an exponential decay of the curvature away 
from the throat into a cone (or wedge) inlet and 
approaching the throat similar to a contour described by a 
circle. 


Step 2: Real Gas Geometric Area Ratio 

The second step is to determine the real gas area ratio 
for the design conditions. A steady, quasi-one- 
dimensional flow field solution that includes all real gas 
effects is needed. This solution can he obtained by either 
solving the steady-state equations directly for a given 
Mach number 0 or by using standard CFD techniques to 
compute the steady-state one-dimensional solution for an 
expanding source flow until the desired exit condition is 
obtained. 

Once the one-dimensional solution for the area ratio 
has been determined, the geometric area ratio for the 
nozzle must be computed. The curvature of the nozzle in 
the throat region introduces a contraction effect due to the 
curvature of the sonic line. Compared to a one- 
dimensional calculation, the contraction effect reduces the 
mass flow through the throat. Si veils’ contraction 
coefficient is computed by using an expansion solution 
for the transonic flow region. 16 For an axisymmetric 
flow 


c () = i-Z±l 
1 96S" 


1 + 


8y ■ 


245 


■2 t 75 y- 


■7577 + 3615 


28805 2 


(4) 


The geometric area ratio is calculated by multiplying the 
one-dimensional area ratio and the contraction coefficient. 


Step 3: Throat to Conical Section Design 

The third step is to determine the contour design from 
the throat to the beginning of the conical section using 
Sivells’ MOC design code. The necessary input 
parameters define the centerline Mach number distribution 
for the MOC design. The input parameters are defined by 
the ratio of the radius of curvature to the throat radius, 
expansion angle, ratio of specific heats, and by 
specifying the Mach numbers at points B and C in Fig. 1 . 
To facilitate the final step, we select Si veils’ option for 
including points along the source flow section and choose 
the axial coordinate to be measured from the upstream 
source for the radial flow region. The gamma is selected 
to be the approximate value for the transonic region 
(Mach = 1.0 to 1.3) and the throat radius is specified to 
be unity. For most gases, the temperature change is 
limited in this range, resulting in only minor variations of 
gamma. 


Step 4: Turning Contour Design 

The next step is to design the contour downstream of 
the conical section by using Sivells’ code for the gamma 
at or near the design Mach number. The same inputs to 
Sivells’ code are used except that the gamma is changed 
and the nozzle exit diameter is specified to achieve the 
geometric area ratio obtained in step 2. 

Step 5; Assembling the Inviscid Nozzle Design 

The final step of the inviscid design is to assemble the 
pieces. First the throat-to-end of the conical section (step 
3) is combined with the turning contour design (step 4). 
Usually a portion of the conical section obtained from the 
throat-to-conical section design overlaps the turning 
contour design. The overlap portion is removed, slightly 
reducing the size of the conical flow region obtained in 
step 3. Next the origin of the axial coordinates is 
transformed from the origin of the source to the nozzle 
throat. Now the subsonic-throat region (step 1) can be 
added directly to the rest of the design. The final design 
coordinates are obtained by scaling both coordinates by a 
common factor to obtain either the desired throat or exit 
radius. 

Limits of Design Procedure 

The above procedure assumes that the design 
parameters selected for the nozzle will result in a conical 
(source) flow region. This can usually be accomplished 
by selecting the appropriate expansion angle. If the first 
choice of expansion angle results in a nozzle design 
without a conical section, the expansion angle can be 
lowered until a reasonable design is obtained. Typically, 
this problem occurs only in low supersonic nozzle 
designs. This procedure also assumes that the gamma 
will increase (or stay constant) as the flow expands. 

Centerline Mach Number Distribution for Use with CFD- 
Based Design Codes or Real Gas MOC Design Codes 
The above procedure was illustrated with the 
assumption that Sivells’ computer code was being used to 
design the different supersonic parts of the nozzle 
contour. There are a few applications where a real gas 
MOC design code or a CFD-based design procedure must 
be used because of either rapidly varying gamma effects 
in the turning region or the BL size. A centerline Mach 
number distribution needs to be prescribed either as a 
boundary condition or as part of the objective function. 
A centerline Mach number distribution for use in the 
design can be determined by using a similar approach as 
Sivells. 

The advantage of Sivells' centerline distribution is that 
it smoothly matches the different flow regions, 
eliminating discontinuities in contour curvature. In 
Sivells’ method, the centerline distribution consisted of 
four sections: (1) a fourth-order polynomial connecting 
the transonic section and the source flow regime, (2) the 
source flow region, (3) a fifth-order polynomial 
connecting the source flow region and the uniform flow 
regime, and (4) the uniform flow region at the design 
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Mach number. The coefficients of the polynomials are 
determined by matching the Mach number or axial 
velocity and its derivatives, with respect to no/zle length, 
at the beginning and end of each section. Sivells derived 
for the ideal gas equations the Mach number derivatives 
with respect to the area ratio for the source flow r region 
and the nozzle sonic line. For a real gas, the coefficients 
of the polynomials and the source flow region need to be 
computed with a real gas model. 

The necessary data for the source flow' and the 
polynomial coefficients can easily be computed with two 
different CFD flow field solutions. The first CFD 
solution is for a two-dimensional (or axisymmetric) flow 
field in the subsonic-transonic region of the nozzle. The 
subsonic contour is defined by Eq. (1) and the supersonic 
side can be described using any expansion that preserves 
continuity at the throat wall in radius, slope, and 
curvature. Because the inviscid supersonic flow field is 
hyperbolic, the supersonic flow field does not affect the 
solution at the sonic line. The throat contraction 
coefficient and the derivatives of Mach number and axial 
velocity at the throat can then easily be computed from the 
flow field solution. This part of the procedure is already 
being used in the design method described by Korte et 
al. n The second solution needed is for a quasi-one- 
dimensional flow field to obtain the Mach number and its 
derivatives with respect to nozzle area ratio. This flow 
field solution is also needed to define the area ratio of the 
nozzle, as described in the section above. A future paper 
will present results for an optimized nozzle design that 
uses CFD-based techniques with an objective function for 
the centerline Mach number as defined in this paragraph. 

Application and Accuracy of the Design Procedure 

An example of the nozzle design process will be 
given in this section, along with computed inviscid CFD 
flow fields to validate the design process. Two nozzles 
will be designed, analyzed, and compared for flow 
quality. The nozzles will be designed with the same 
parameters except that the ratios of specific heats will be 
different to account for different gas assumptions. 

One nozzle will be designed for an ideal gas by using 
Sivells’ design code with a gamma computed to give the 
design Mach number at the same area ratio as the real gas 
design. The other case wall use the real gas procedure for 
a calorically and thermally imperfect gas. The same 
parameters will be used for defining the centerline Mach 
number distribution, except for the ratio of specific heats. 
The flow quality of the nozzle designs will be evaluated 
based on the CFD inviscid simulations. 

Nozzle Design 

To demonstrate the procedure, we will design a 
nozzle for an exit Mach number of 14, with stagnation 
conditions that have significant real gas effects. These 
specifications are made based on the availability of a 
validated real gas equation of state for high pressure and 
temperature nitrogen from previous work. 1 * 


For the flow conditions given in Table 1, the one- 
dimensional real gas area ratio is computed to be 
2396.87. If only caloric imperfections are considered, 
then the area ratio is computed to be 3232.44. The ratio 
of specific heats (y) for an ideal gas needed to obtain 
Mach 14 and an area ratio of 2396.87 is 1.4073. If the 
test gas is air and it is treated as an ideal gas, the area ratio 
for Mach 14 with y = 1 .4 is 2686. The nonideal effects of 
the high pressure (on the thermal equation of state) and 
the high temperature (on the caloric equation of state) 
offset each other so that the real gas area ratio is only 
10.8% lower than the ideal gas area ratio for y= 1 .4. 


Table 1. Nozzle design conditions and parameters 


Total pressure 

20,000 psi 

Total temperature 

3,000°R 

Gas composition 

Nitrogen 

Stagnation compressibility factor 

1.31205 

Exit Mach number 

14.0 

Subsonic approach angle 

45° 

Ratio radius-of-curvature/throat radius 

3.0 

Supersonic (source flow) angle 

12.0° 

Throat radius 

1.0 

Geometry 

Axisymmetric 

Table 2. Area ratios and ratio of specific heats 

Case Gas model Mach y 

Area ratio 


1 

Real gas 

1.0 

1.3180 


<( 


1.3 

1.3234 


u 

it 

14.0 

1 .4000 

2396.87* 

2 

Ideal gas 

14.0 

1 .4073 

2396.87 


*For conditions given in Table 1 . 


Two nozzles were designed for the specifications 
given in Tables I and 2. The input for Sivells’ code for 
each case is given in the appendix Table A 1 . The real gas 
area ratio is computed with a real gas equation of state for 
nitrogen. The ideal gas nozzle design is only for 
comparing flow quality with that of the real gas design. 
The subsonic contour is identical for both nozzle designs 
(appendix, Table A2). The supersonic contours for both 
nozzles arc given in appendix Table A3. 

CFD Analysis of Nozzle Design 

A CFD flow field solution for each nozzle was 
computed by using the methods given in Refs. 1 i and 
18. Both the designs and the CFD analyses used the 
same equation of state. The throat region was computed 
by solving the unsteady Euler equations. The throat 
numerical grid was 81 x 75 points in the axial and normal 
directions, respectively. The centerline solution for the 
throat started at approximately Mach 0.2 and ended at 
approximately Mach 1.8. The downstream portion of the 
nozzle was computed by solving the steady Euler 
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equations with a space marching method and was 
computed in approximately 1000 marching steps. 

The computed centerline Mach number distributions 
for the whole nozzle in each case are compared in Fig. 3. 



Fig. 3 Comparison of computed centerline Mach 
number distributions for real and ideal gas nozzle 
designs. 

At this scale there are no differences in the two designs, 
and both appear to generate the ideal uniform core region 
at Mach 14. Focusing in on the uniform core region in 
Fig. 4 shows slight Mach number variations on the order 
of 0.02 Mach for each case. 



Fig. 4 Comparison of computed centerline Mach 
number distributions for real and ideal gas nozzle 
designs in uniform core region. 



Fig. 5 Comparison of exit Mach number profile for real 
and ideal gas nozzle designs. 


Real Gas Design (Case 1) 



Fig. 6 Comparison of exit flow angle profile for real and 
ideal gas nozzle designs. 

The exit Mach number profiles are given in Fig. 5, and 
the angle profiles are given in Fig. 6. The flow is 
essentially uniform at both nozzle exits because the 
magnitude of the flow angles is so small. The Mach 
number profile at the exit varies approximately 0.02 Mach 
across the nozzle height. The computed exit Mach 
numbers for both nozzles were slightly higher than the 
Mach 14 design point: Mach 14.04 for the ideal gas 
design and Mach 14.02 for the real gas design. This 
minor variation from the design point is most likely due 
to the differences in the transonic expansion solution for 
the nozzle throat used in the design. Mach contours for 
the real gas nozzle design shown in Fig. 7 show only a 
0.01 Mach number variation in the uniform core region. 



Fig. 7 Computed Mach contours for a real gas nozzle 
design. 


In both designs, minor variations from the Mach 14 
design point are observed. However, the variations arc 
similar and are most likely characteristic of Sivells' 
design process, not of the real gas design procedure 
presented above. In fact, the computed solution for the 
real gas nozzle has flow quality almost identical to that of 
the computed ideal gas nozzle design. Variations of this 
real gas design procedure have been used to design 
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nozzles at Langley and elsewhere with excellent 
results. 1521 

Summary 

A simple, straightforward method for the inviscid 
design of a real gas nozzle for hypersonic facilities has 
been proposed. The design is demonstrated for real gas 
conditions with significant calorically and thermally 
imperfect effects. An ideal gas nozzle is also designed 
for the same area ratio and Mach number for comparison 
of flow quality. CFD flow field solutions are computed 
for each nozzle, and the flow field qualities in terms of 
Mach number variation and flow' angle arc compared. 
The nozzle designs have nearly identical computed flow 
quality. This method is easy to use and is extremely 
accurate when the real gas effects occur before the 
beginning of the turning contour design. 
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Table A 1 . Input data for Si veils’ design code* 

Case 1: Real gas design, upstream section-supersonic. 

12345678901234567890123456789012345678901234567890123456789012345678901234567890 
R_Gas-l 0 


1.323 

1716.563 

1.000 


.896 

0.227E 

:-07 

.000 

.000 

1000.000 

12.000 

3.000 

.000 

11 

.200 

14. 

000 

.00000 

- 1.000 

10.000 

41 21 

0 10 

1 41 

49 

61 

32 

0 

-1 -1 

10 0 

-21 13 

Case 1 : Real 

gas design, 

downstream i 

section-supersonic . 



R Ga s — 1 

0 









1.400 

1716.563 

1.000 


.896 

0.227E 

-07 

.000 

.000 

1000.000 

12.000’ 

3.000 ’ 

.000 

11 

.200 

14. 

000 

-48.91131 

- 1.000 

10.000 

41 21 

0 10 

-1 41 

49 

61 

32 

0 

-1 -1 

10 0 

-21 13 

Case 2: Ideal gas design, supersonic 

contour . 





I -Gas -2 

0 









1.4073 

1716.563 

1.000 


.896 

0.227E 

-07 

.000 

.000 

1000.000 

12.000 

3.000 

.000 

11 

.200 

14. 

000 

.00000 

.000 

10.000 

41 21 

0 10 

0 41 

49 

61 

32 

0 

-1 0 

10 0 

-21 13 

*input data for 

running inviscid design only. 








Table A2. Subsonic contour 


X 


y«' 

- 4.9462 

4.2087 

- 1.0000 

- 4.8472 

4 . 1098 

- 0.9996 

- 4.7483 

4.0109 

- 0.9984 

- 4.6494 

3.9123 

- 0.9963 

- 4.5505 

3.8138 

- 0.9934 

- 4.4515 

3.7158 

- 0.9897 

- 4.3526 

3.6181 

- 0.9851 

- 4.2537 

3.5209 

- 0.9796 

- 4.1548 

3.4243 

- 0.9732 

- 4.0559 

3.3284 

- 0.9659 

- 3.9569 

3.2332 

- 0.9578 

- 3.8580 

3.1389 

- 0.9487 

- 3.7591 

3.0455 

- 0.9387 

- 3.6602 

2 . 9532 

- 0.9278 

- 3.5612 

2.8620 

- 0.9160 

- 3.4623 

2.7720 

- 0.9033 

- 3.3634 

2.6833 

- 0.8897 

- 3.2645 

2.5960 

- 0.8752 

- 3.1655 

2.5102 

- 0.8598 

- 3.0666 

2.4259 

- 0.8435 

- 2.9677 

2.3433 

- 0.8263 

- 2.8688 

2.2625 

- 0.8082 

- 2.7699 

2 . 1835 

- 0.7893 

- 2.6709 

2.1064 

- 0.7695 

- 2.5720 

2.0313 

- 0.7489 

- 2.4731 

1.9582 

- 0.7275 

- 2.3742 

1.8873 

- 0.7053 

- 2.2752 

1.8187 

- 0.6823 

- 2.1763 

1.7524 

- 0.6585 

- 2.0774 

1.6885 

- 0.6340 

- 1.9785 

1.6270 

- 0.6088 

- 1.8795 

1.5680 

- 0 . 5829 

- 1.7806 

1.5117 

- 0.5563 

- 1.6817 

1.4580 

- 0.5291 

- 1.5828 

1.4070 

- 0.5013 

- 1.4838 

1.3588 

- 0.4729 

- 1.3849 

1.3135 

- 0.4439 

- 1.2860 

1.2710 

- 0.4144 

- 1.1871 

1.2315 

- 0.3845 

- 1.0882 

1.1950 

- 0.3540 

- 0.9892 

1.1615 

- 0.3232 

- 0.8903 

1.1310 

- 0.2920 

- 0.7914 

1.1037 

- 0.2604 

- 0.6925 

1.0795 

- 0.2286 

- 0.5935 

1.0585 

- 0.1964 

- 0.4946 

1.0407 

- 0.1640 

- 0.3957 

1.0261 

- 0.1315 

- 0.2968 

1.0147 

- 0.0987 

- 0.1978 

1.0065 

- 0.0659 

- 0.0989 

1.0016 

- 0.0330 

0.0000 

1 . 0000 

0.0000 
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